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We study the influence of the phonon environment on 
the electron dynamics in a doped quantum dot molecule. 
A non-perturbative quantum kinetic theory based on cor- 
relation expansion is used in order to describe both diag- 
onal and off-diagonal electron-phonon couplings repre- 
senting real and virtual processes with relevant acoustic 
phonons. 



1 Introduction Quantum dots (QDs) are one of the 

promising candidates for building a feasible quantum com- 
puter. In particular, it has been proposed to use two coupled 
quantum dots called quantum dot molecules (QDMs) for 
various schemes of quantum computation, where the two 
ground states of a single confined electron as well as sin- 
glet and triplet states of doubly doped structures can be 
used as the logical qubit states. The electrical G) and op- 
tical control lEE] °f spin- and charge-based qubits have 
already been demonstrated. However, in such solid state 
systems, the phonon-assisted relaxation (5][6] can strongly 
affect the coherent control. The energy difference between 
the two electron ground states, which is typically of order 
of a few meV, can lead to pure dephasing processes JVj[Sj 
9 1 and to electron tunneling I10II1 1II12L 

In this paper, we present the full quantum kinetic 
description of the phonon-mediated relaxation in doped 
quantum dot molecules including non-Markovian effects. 
We show that the coupling to the phonon reservoir in quan- 
tum dot molecules can lead to a fast electron tunneling on 
a picosecond timescale, which strongly affects the co- 
herent electron evolution. We employ the non-Markovian 
correlation expansion technique H13U14I including up to 
three-particle correlations. Due to space constrains, we 



We show that the relaxation is dominated by phonon- 
assisted electron tunneling between constituent quantum 
dots and occurs on a picosecond time scale. The depen- 
dence of the time evolution of the quantum dot occupa- 
tion probabilities on the energy mismatch between the 
quantum dots is studied in detail. 
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present here equations of motions with up to two-particle 
correlations. We include diagonal and off-diagonal cou- 
plings to the phonon reservoir representing virtual and real 
phonon-assisted processes, respectively. The dependence 
of the time evolution of the quantum dot occupation prob- 
abilities on the energy mismatch between the constituent 
quantum dots is studied in detail. 

2 Model system We consider a system consisting of 
a single quantum dot molecule doped with one electron. 
The relaxation between the two energetically lowest states 
of the electron 1 1) in the left and 2} in the right quantum 
dot is considered (see Fig. 1 with schematic plot of the en- 
ergy levels in the quantum dot molecule). The free Hamil- 
tonian of the electron reads 

H c = e (|2)(2|-|l)(l|) + r(|l)(2| + |2)(l|), (1) 

where r is the direct tunneling coupling element between 
the quantum dots and Ae = 2e is the energy difference 
between the ground states in both quantum dots. 

The Hamiltonian describing the free evolution of 
phonons is 



H, 



P h 



J2 htj 



bi .bu 



(2) 
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Figure 1 A schematic plot of a quantum dot molecule with 
two electron states |1) and 2) with the energy difference 
of Ae — 2e. The dashed and curled arrows depict electron 
tunneling with simultaneous phonon emission. 



with 6^ and 6k, s being phonon creation and annihila- 
tion operators, respectively. cjk,s denotes phonon frequen- 
cies with phonon wave vector k, where different phonon 
branches are labeled by s (longitudinal I and two transverse 
tl and t2). 

The interaction between the electron and phonon reser- 
voir reads 



H; nt = 



l mt — / 

i,j=l,2 k,s 



(3) 

The relevant coupling elements for the acoustic phonons 
coupled via both piezoelectric coupling and deformation 
potential read: 



2pc s k £q£ t 



(4) 



with Gaussian electron wave functions 

^ (r) ^^p cxp (-w 



D e I d 3 r^(r)e lk ' r ^(r) (5) 



(6) 



We have numerically simulated the self-assembled GaAs 
quantum dots with the following parameters: l\ = 4 nm, 
1-2 = 4.1 nm, and d = 6 nm, where ij is the electron 
wave function size of the i-th quantum dot, and the dis- 
tance between the dots is d. Here, p = 5360 kg/m 3 is 
the crystal density, V is the normalization volume of the 
phonon modes, c s is the speed of sound (longitudinal q = 
5150 m/s or transverse c t = 2800 m/s, depending on the 
phonon branch), dp = 0.16 C/m 2 is the piezoelectric con- 
stant, D e = —8 eV is the deformation potential constant 
for the electrons and e r = 13.2 is the static dielectric con- 
stant. The function M s (k) in the piezoelectric coupling el- 
ement reads: 

Mi{k) = \ sin(0) sin(20) sin(2^), (7) 



M a {k) = -sin(26»)cos(2v?), 

M t2 {k) = sin(6»)[3cos 2 (^) - 1] sin(2</>). 



(8) 
(9) 



We study the time evolution of the electron confined 
in a quantum dot molecule in the presence of the coupling 
to the phonon reservoir within the density matrix theory. 
We employ the second order correlation expansion method 
01311141 . where it is assumed that correlations involving an 
increasing number of particles are of decreasing impor- 
tance. This non-perturbative technique covers the memory 
effects in the non-Markovian regime. 

With the help of Heisenberg equation we derive the 
equations of motion for the quantities of interest. The equa- 
tion for the electron occupation / = (|2)(2|) reads: 



-r Im(p) 



k 

E^( k ) [' 



Sk 



„(+)* 



,(+) 



2iB k Im(p) 



(10) 



It couples to the coherence p = (|1)(2|) as well as to 
the phonon-assisted coherences s q = (|l)(2|6 q ) cor '' and 

«q +) = (|l)(2|6 q ) corr and to phonon amplitudes B k = 
(6k). The factorization scheme has been used above with 



(|l)(2|6 k ) = (|l)(2|)(&k) + (|l)(2|&k) c 



(11) 



where the quantities have been decomposed into all pos- 
sible lower-order factorizations. The occupation is directly 
affected by phonon reservoir in the case of the off-diagonal 
coupling ~ 5i2 (k) and indirectly by the coherence p, 
which evolves in the following way: 



P 



k 



^(1-2/) 
5l2 (k) [2t k + (1 - 2/)5 k ] 



(12) 



~5>ia(k) [2t* k + (1 - 2f)B* k ] 



k 



[322 (k) -gii (k)] (sk+pS k ) 



3n(k)] 



,(+) 



Here, a direct impact from phonons is present not only via 
real transitions but also for the pure dephasing described 
by the diagonal coupling elements <?n(k) and 322 (k). If 
these coupling elements are equal, the two pure dephasing 
channels cancel each other out . In this case, the informa- 
tion cannot leak out, since the phonons is not delivering 
any which way information and no decoherence channel 
exists. 

The evolution of phonons reveals the diagonal and off- 
diagonal couplings to the density and coherence, respec- 
tively, 



(13) 



-[<? 2 * 2 (q)-<?n(q)](i-/)- 
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The next step is to derive the equations of motion for 
all the phonon-assisted correlation quantities. 



s' q = -j- 2rt <i - ft [-922 (q) - 9*11 (q)] fp 

i 



(14) 



,(+) 



-(2e + 7iw q ) Sq - -fff 2 (q)[l - / - 2Re(p)p] 



J$>i 2 (k) [(1 - 2/)n kq ) + 2t q B k 



jE-9i*2( k ) [(l-2/)n kq + 2t q sa 

k 

-££ [922 W - 911 (k)] (pn^ } + s q 5 k 



4 £ [5 2 * 2 (k) - 5u(k)] (pn kq + s qJ B k ) , 



~2rt;-i fll2 (q)[2Re(p)jj-/] 
--(2e - Tio; q )s q +) - -|>22(q) - ffii(q)]/p 
-l^ ffl2 (k) [(1 - 2/)n qk + 2t*B k ] 



(15) 



£$>* 2 (k) (l-2/)n k ; ) * + 2*;S k 



4 E t92 2 (k) - ^(k)] (pn£* + Sq +)B k 



The dynamics of the phonon-assisted density, i q 
{|l)(l|fe q ) corr is described by 



-i(w k + w q )n qk ) - (q) 



Sk+ 4 +) * 



(18) 



"9^2 (k) 



„(+)* 



[«? 2 * 2 (k)- ffl \(k)]i c 



+^ [922 (q) -9u(q)]*k. 



In addition, the two-particle correlations couple to three 
particle correlations, e.g. (|l)(2|6 q s b 



k,s' 



. and these 

couple to up to four particle correlations, etc. Thus, in or- 
der to get a closed set of equations, one needs to truncate 
the hierarchy by neglecting higher order correlations. In 
the present paper, up to three particle correlations were in- 
cluded. 

3 Phonon-mediated tunneling Initially, one elec- 
tron is injected into the quantum dot with higher energy, 
/ = (|2)(2|) = 1. The time of the injection is assumed to 
be much shorter than the response of the phonon reservoir, 
thus at the initial time we can set all correlations to zero. 

The time evolution of the occupation probability of the 
quantum dot with higher energy in the presence of the cou- 
pling to the phonon reservoir is shown in Fig. [2] for differ- 
ent values of the energy difference Ae between the two 
electron ground states. The case of the ideal evolution, 
where the electron-phonon coupling is absent, is shown 
in the inset of Fig. [2] The electron tunnels coherently be- 
tween the two quantum dots with a period of these oscilla- 
tions that is determined by the energy difference between 
the ground states Ae = 2e and the tunneling coupling r 
and scales like ~ r/(Ae). The maximum occupation also 
depends on these parameters. If the energy difference is 
twice the value of the tunneling coupling, the electron tun- 
nels between the quantum dots with maximum occupation 
of \ of the quantum dot with lower energy. For the larger 



tq = -*Wqtq - T9*M [P ~ 2Re(p)(l - /)] 



(16) 



-r\ s -*(+>* 



+ ft[92 2 (q)-9n(q)]/(l-/) 



■~ J2 9i 2 (k) [(«, - Sq +) *) B k + n^2ilm(p) 



~ £ 9i2(k) [(«„ - Bl + n kq 2zlm(p) 



All these phonon-assisted two-particle correlations couple 
to the phonon-phonon correlations, phonon density n qk = 



(bqb]s) corr and phonon coherence n kq = (& k & q ) corr : 



« q k = -*(wk - w q )ri qk + -c/i 2 (q) 



, (+)* 

sk + 4 



(17) 



9i* 2 (k) 



{g; 2 (k) - g* n (k)]t* c 



[522 (q) - 511 (q)] *k, 




60 
time (ps) 



Figure 2 The occupation probability of a quantum dot with 
higher energy as a function of time for different values 
of the energy difference Ae between two electron ground 
states for a fixed tunneling coupling r = 0.1 meV and at 
fixed temperature T = 4 K. 
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energy differences, the tunneling coupling is too weak to 
drive the electron into the adjacent quantum dot and the 
electron stays in the initial quantum dot. 

The dynamics strongly changes if the coupling be- 
tween the electron and the phonon reservoir is present (see 
Fig.[2]i, which results in the electron tunneling to the neigh- 
boring quantum dot on a short picosecond time scale. This 
phonon-assisted tunneling is only properly modeled if the 
off-diagonal electron-phonon coupling term in the Hamil- 
tonian is included. Both coupling via deformation poten- 
tial and piezoelectric effect are important and contribute 
to the relaxation. In the case of the energy difference of 
Ae = 0.2 meV, the thermalization process takes about 
50 ps and finishes with a delocalized state of the electron 
with a probability of ~ 0.7 for being in the left quantum 
dot and ~ 0.3 in the second dot. For larger energy mis- 
matches of 1 and 2 meV, the relaxation takes about 30 and 
40 picoseconds, respectively, and in the final state, the elec- 
tron is localized in the quantum dot with lower energy. In 
these cases, the probability of phonon emission or/and ab- 
sorption processes is the highest since the energies of the 
relevant acoustic phonons lie in this parameter regime. If 
we increase the energy difference between the dots, the 
probability of phonon-mediated relaxation will decrease 
and the tunneling is slower and less efficient. The presented 
results were calculated for the low temperature of T = 4 K. 
The phonon-assisted tunneling will get stronger and faster 
at higher temperatures and will also affect the final state, 
since the ratio of probabilities of phonon emission and ab- 
sorption changes with temperature. 

4 Conclusion A full description of the electron dy- 
namics in the presence of the electron-phonon coupling 
in a quantum dot molecule doped with a single electron 
has been presented. The two couplings to the relevant 
acoustic phonons via deformation potential and piezoelec- 
tric coupling with diagonal and off-diagonal interactions 
have been taken into account. It has been shown that the 
phonon-mediated relaxation is a fast process on a picosec- 
ond timescale strongly modifying the coherent evolution 
of the electron. It is dominated by off-diagonal electron- 
phonon coupling. We analyzed the dependence of the 
relaxation on the energy difference between the two quan- 
tum dots and indicated the vales of parameters when the 
phonon-mediated tunneling is most efficient. 
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